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Kubo formula is used to get the d.c conductance of a statistical ensemble of two dimensional clusters 
of the square lattice in the presence of random magnetic fluxes. Fluxes traversing lattice plaquettes 
are distributed uniformly between — and with $o the flux quantum. The localization 
length is obtained from the exponential decay of the averaged conductance as a function of the 
cluster side. Standard results are recovered when this numerical approach is applied to Anderson 
model of diagonal disorder. The localization length of the complex non-diagonal model of disorder 
remains well below £ = 10 4 in the main part of the band in spite of its exponential increase near 
the band edges. 

PACS: 72.15.Rn, 71.55.Jv, 73.40.Hm, 74.72.-h 
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I. INTRODUCTION 

The two-dimensional (2D) motion of a single particle 
in a random magnetic field is attracting strong inter- 
est since the model is related both to. the slave-boson 
description of high T c superconductorsEI and the Chern- 
Simons theory of the half-filled Landau level.0 The ex- 
istence of recent experiments measuring transport prop- 
erties in a static random magnetic field adds consider- 
able interest to this subject.13 In principle, all states of 
a 2D system are exponentiallyJocalized according to the 
scaling theory of localizations The presence of a mag- 
netic field that destroys time-reversal symmetry makes 
a stand against localization. Nevertheless, perturbative 
renormalization gwup calculations show that all states 
remain localized Bu 

While theory points towards localization, numerical 
work has not led to a consensus of opinion.EI Loosely 
speaking, states look extended in finite samples as soon 
as the edge of the band is left and a definite answer is 
quite difficult. Therefore, the most extended conclusion 
is that a mobility edge separates localized states from 
extended or critical states. 

The model of 2D random fluxes has been previously 
studied from a statistical point of view. The analysis of 
several statistical magnitudes (nearest-neighbor spacing 
distribution, the probability distribution of wavefunction 
weights, etc.) shows a systematic flow from the Gaus- 
sian Unitary Ensemble (GUE) statistics at small system 
sizes to Poisson statistics at larger sizesfl If GUE and 
Poisson statistics are the only possibilities in the macro- 
scopic limit, localization over the whole spectrum follows 
from the referred numerical results. If, on the contrary, 
more fixed points are possible, then the flow could even- 
tually start near the metallic universality class and stop 
at certain new class characterized probably by some frac- 
tal properties (for example, some power law increase of 
the number of sites visited by an eigenstate). 

Both the theoretical approaches to the metal-insulator 
transition and the numerical tools used for the study 



of localization have a long history.H'tTO The standard 
numerical method has been really successful describing 
the critical behaviot-pf the 3D metal-insulator transition 
of Anderson model. O Nevertheless, the numerical proce- 
dure is somewhat artificial. Let me remind that firstly, 
localization lengths of ID systems (strips or bars) of in- 
creasing transverse size are calculated and, secondly, ex- 
trapolation of the data gives the localization length of 
infinite 2D or 3D systems. This introduces a fictitious 
symmetry breaking of the 2D or 3D system since a partic- 
ular direction is preferentially treated. Moreover, being 
all ID states exponentially localized, they become even- 
tually extended after doing finite-size scaling analysis of 
the data. This last part of the calculation introduces 
some uncertainties in the procedure. 

In this work, I return to the original idea of meas uring . 
the conductance of finite samples on the computerli3'E2l 
Today's use of the older scheme has some advantages: 
mesoscopic physics has led to a more efficient use of Kubo 
formulali-j, computers are much more powerful, there is 
a large experience in the numerical calculation of local- 
ization lengths, etc.. Based on the statistical analysis 
of conductance data obtained for samples of increasing 
size, the direct evaluation of the localization length is 
possible. I find this approach really appealing for two 
reasons: i) the exact static conductance is obtained via 
Kubo formula (neither a small imaginary part should be 
added to the energy nor other numerical trick should be 
used to reach the static limit of the conductance of a fi- 
nite sample), and ii) the scaling behavior of a measurable 
magnitude (conductance) is directly analyzed. 

The paper is organized as follows. The model that is 
numerically studied in the paper is defined in Section II, 
some remarks about the use of the Kubo formula within 
the present context are given in Section III, the size de- 
pendence of the conductance is discussed in Section IV, 
results for the Anderson model and the model with ran- 
dom fluxes are presented in Section V, and finally, Section 
VI ends the paper with some remarks. 
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II. LATTICE MODEL 

The Hamiltonian describing random magnetic fluxes 
on a L x L cluster of the square lattice is : 

L 2 

# = 5>c+q- i£e 2 ^»'4c, , (1) 
;=i <w> 

where cj creates an electron on site I, I and V are nearest- 
neighbor sites, e; are random diagonal energies chosen 
with equal probability from the interval [— W/2,W/2], 
and —t is the hopping energy (hereafter t = 1 is cho- 
sen). The flux through a given loop S on the lattice is 
$s = X)<; 2'>es c t )u ' measured in units of flux quantum 
$ = — • Link fields satisfy 4>w — ~4>i'i- Either link 
fields (Meissner phase) or magnetic fluxes (Debye phase) 
can be randomly selected.E£l Here, the second possibil- 
ity has been chosen: uncorrelated fluxes are randomly 
selected with equal probability from the interval [— ^, ^] 
while W = 0. Notice that we are considering the model 
with pure non-diagonal disorder that shows maximum 
magnetic disorder. This model will be called RMF (Ran- 
dom Magnetic Fluxes or Random Magnetic Field) model 
hereafter. 

Some additional calculations have been done for a real 
version of Hamiltonian (Q). As before, diagonal disorder 
is absent (W = 0) but now link field values are restricted 
to and 1/2 that are randomly chosen with equal proba- 
bility. The model describes a square lattice in which the 
hopping integral is constant in absolute value but takes 
random signs. This model will be shortly called RHS 
(Random Hopping Signs) model. Finally, a pure Ander- 
son model of real diagonal disorder has also been studied 
in order to compare with results in the literature. In this 
case, W ^ and all link fields vanish. 

III. D.C. CONDUCTANCE VIA KUBO FORMULA 

According to Kubo's linear resppase theory, the static 
electrical conductivity is given by:E3 

G = a xx (0) = -(^-)2Tr[{hv x )Img(E)(hv x )lmg{E)\ , 
h 

(2) 

where the velocity (current) operator v x is related to the 
position x operator by the equation of motion: 

ihv x = [H, x] (3) 

and ImC? (E) is defined in terms of advanced and retarded 
Green functions as: 

hng(E) = hg+(E)-g-(E!)] . (4) 

The conductance is calculated for a system consisting 
of a L x L cluster of the square lattice connected to two 



ideal semiinfinite leads of width L. The exact form of 
the electrical field does not matter in a linear response 
one-electron theory. Therefore, an abrupt potential drop 
at one of the two cluster sides provides the simplest nu- 
merical use of the Kubo formula. In this case, operator 
v x has finite matrix elements on only two adjacent layers 
(see Eq.(g)) and Gr een functions are just needed for this 
considerably restricted subset of sites owing to the trace 
appearing in Eq. (^) . Since Green functions are given by: 

[E-H-£ T (E)-± 1 (E)]g(E)=I , (5) 

where S r (-E) and l£\(E) are the selfenergies introduced 
by the semiinfinite right and left leads, respectivelyli-3, 
the evaluation of the desired elements of ImO(E) is ef- 
ficiently achieved applying an LU decomposition to the 
band matrix 

<l\E-H-t T (E)-ti(E)\l' > . (6) 

Selfenergies are troubleless calculated for any energy E 
using a recurrent algorithm. Retarded or advanced Green 
functions (g + (E) or Q~(E)) are obtained using the cor- 
responding selfenergies in Eq.(|5|). Since lead selfenergies 
add an imaginary part to some diagonal elements of the 
Hamiltonian, there is no need to add a small imaginary 
part to the energy E before solving the set of linear al- 
gebraic equations. 

The discussion of a technical point is here in or- 
der. It has been proved that Kubo formalism leads to 
the more familiar Landauer-Buttiker expression of the 
conductance as a sum of transmission coefficients over 
channels .Ell Usually, transmission is obtained within a 
transfer matrix technique in which the selfenergy pro- 
duced by the lead at the left is iterated through the dis- 
ordered part of the sample and matched at the right side 
to the selfenergy coming from the right lead. In this cal- 
culation, matrices of dimension L are inverted L times. 
It could seem that this method is computationally more 
appropriate because it deals with smaller matrices (the 
elements of the Green function within the sample are 
given by the inverse of matrix (|6|) which has dimension 
L 2 ). Nevertheless, since matrix (||) is actually a band 
matrix of bandwidth 2L + I and only a few elements of 
the Green function are necessary for the evaluation of the 
conductance via Kubo formula, standard numerical pro- 
cedures (LU decomposition of a band matrix followed 
by both a forward substitution and a backsubstitution 
for each needed Green function element) are quite ad- 
vantageous. Moreover, canned subroutine packages cus- 
tomized for particular computers can be used. Thus, we 
end with a very simple computer code written in a high 
level language that can be advantageously run on any 
given computer. Typically, 127 seconds are needed for 
the computation of the conductance of a 100 x 100 sam- 
ple on an alpha DECstation 3000. 
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IV. SCALING OF THE CONDUCTANCE 

Conductance is measured for a statistical ensemble of 
L x L clusters with random disorder. This provides a set 
{Gi} of N values of the static conductance that allows 
the evaluation of both the average value of conductance 



1 N 

»=i 



(7) 



and its variance 



1 N - 

Var(G 1 ...G w ) = ^- T ^(G l -G) 2 

i=l 



(8) 



The error bar of G is given by its standard deviation 



< G 2 > -< G > 2 



which is well approximated by 



<7(G) 



Var(Gi-Giv) 
N 



(9) 



(10) 



as soon as iV is large enough. Since conductance shows 
typical universal fluctuations of order 1 (in units of e 2 /h), 
the error of the measured average conductance is ~ N~?. 

This way of characterizing the whole distribution by 
just two values (G and its error bar) works fine for mag- 
nitudes showing an almost Gaussian statistics. Never- 
theless, the distribution of conductances clearly deviates 
from normal when the regime of exponential localization 
is being approached. In this regime both theory and 
numerics suggest that the distribution of the logarithm 
of the conductance is the appropriate statistical variable 
showing a normal distribution around its average value. 
Therefore, the usual practice is the analysis of logG 
whenever the distribution of G deviates significantly from 
normal. I am not going to follow this standard way of 
analyzing the conductance for several reasons: i) In prin- 
ciple, since G is measured, it seems preferable to work di- 
rectly with it. Moreover, G is both bounded from below 
(0) and from above (L, the number of channels). While 
conductance will vanish in the localization regime, its log- 
arithm will diverge. Consequently, the variance of G is 
bounded but the variance of logG is not. Computation- 
ally, this means that logarithm averages show more noise 
in finite sampling, ii) The statistical characterization of 
the conductance should work from the metallic regime to 
the localized one since regime can change as a function 
of the system size. The switch from G to log(G) is some- 
what ambiguous, iii) The median plays quite good the 
role of a central value of some statistical magnitude even 
when the distribution of this magnitude is not Gaussian. 
For example, I have checked that the conductance me- 
dian G me d gives exactly the same exponentially decaying 
behavior as log(G) in the localization regime. Actually, 



the median will work for any variable change that pre- 
serves ordering among data, i.e., if /(Gi) < f(Gj) when- 
ever Gi < Gj . This property is quite general and makes 
the use of the median quite robust, iv) Unfortunately, 
the median has two statistical disadvantages. First, even 
for a normal distribution of the original variable, its vari- 
ance for a finite sampling is larger that the corresponding 
variance of the mean value.Efl Second, the estimation of 
its error bar is not as straightforward as the calculation 
of the standard deviation of the mean (Eq.|l0|). Actu- 
ally, one is forced to divide the database into subsets, 
obtain several median estimates and calculate their vari- 
ance. The result is a bad estimation of the uncertainty of 
the median, v) In conclusion, I will use the conductance 
median without error bar as the magnitude giving the 
central value of the conductance distribution whenever 
a transition from the metallic regime to the localization 
regime takes place somewhere within the computation. 
On the other hand, the use of the mean value and its 
error bar will be preferred when the conductance distri- 
bution remains normal. From here on, either G or G me d 
will be referred as conductance G and will be given in 
units of the conductance quantum e 2 /h. 

Localization is an asymptotic property of the conduc- 
tance as a function of the linear cluster size L. In order 
to correctly describe this behavior using measurements 
done on finite samples, finite size corrections should be 
allowed. Notice that the conductance of a 2D square clus- 
ter of side L increases linearly in the absence of disorder; 
in fact, the conductance is an integer value that gives the 
total number of channels that are open at the Fermi en- 
ergy (this number is exactly L at the band center). The 
following scaling models have been used in this work: 

or 

G(L) — (a + bL^ 1 + cL~ 2 + ...) exp( — — ) , (11) 



and 



G(L) = -^-exp(-^). 
a + oL £ 



(12) 



In both cases, the term multiplying the exponentially de- 
caying factor is designed to describe finite size effects. £ 
measures the localization length. The first expression is 
more intuitive and works always fine (with a large enough 
number of parameters) while the second expression de- 
scribes better the metallic regime in which the conduc- 
tance increases linearly with the system size L. All the 
results presented in the next Section allow a good fit with 
just three parameters (a, b, £). The localization length 
does not depend sensibly on the used model. 

Finally, let me comment about the reliability of the 
localization length obtained by the finite scale modeling. 
Generally speaking, parameters a, 6, c, ... are determined 
by conductance values of small cluster sizes while the 
localization length £ is determined by the behavior of 
conductance at larger system sizes. Since conductance 
starts decreasing after a region of initial increase, the es- 
timation of the localization lengths is reliable only if the 
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crossover to the region in which conductance decreases 
can be reached by numerical simulation. According to 
the data presented in this paper, this happens for lo- 
calization lengths that are well below 10 (lengths are 
everywhere given in units of a, the lattice constant). No- 
tice that the existence of localization at large L is built 
in the model (see Eqs.(|il"|-|l2|)): a decaying conductance 
does not necessarily mean an exponentially decaying law. 
In other words, numerical simulation does not prove lo- 
calization but gives a measure of its characteristic length 
£. Nevertheless, after trying a lot of alternative (and 
worser, i.e., with larger \ 2 deviations from data) mod- 
els for the description of the conductance dependence on 
system size, one gets more or less convinced that-expo- 
nential decay at large L values is quite plausibleJij 

V. RESULTS 

Although my primary interest is the determination of 
localization lengths for Hamiltonians describing random 
magnetic fluxes, the paradigmatic Anderson model of di- 
agonal disorder has also been studied in order to check 
the numerical tools used in this work. Conductance val- 
ues are everywhere given in units of e 2 /h. 

A. Anderson model 

Hamiltonian ([!]) without link fields (<f>w — for all I, V 
pairs) describes the Anderson model of diagonal disor- 
der. The width W of the diagonal energies distribution 
measures the degree of disorder. Exponential localiza- 
tion of all wavefunctions is the accepted state of the art 
for this model. Nevertheless, measuring the localization 
length when disorder is weak (W < 2) is numerically 
very difficult. Wavefunctions are clearly localized in the 
band tails but look extended well within the band. This 
fact has originated more than one claim pointing to the 
existence of a metal-insulator transition at some critical 
energy (mobility edge) within the band. 

Following tradition, I have collected conductance data 
at the band center (E = 0). It is assumed that states 
at finite energies are more localized than the band center 
one. Fig. 1 shows the behavior of both the mean and the 
median of the conductances samples as a function of the 
side L of the L x L disordered cluster in which the cur- 
rent is measured. Since disorder is very large (W = 10), 
both measures of a central value show a clear exponen- 
tial decay. Nevertheless, there are also important dif- 
ferences: localization lengths do not coincide and mean 
values show more dispersion, i.e., they do not follow a 
straight line as good as the median does. Looking closer 
to the data, the reason for the disagreement is obvious: 
conductance is very small (G <C 1), its distribution is 
not normal and, consequently, the mean is dominated by 
outliers conductances of a relatively large value. On the 



other hand, the logarithm of the conductance seems to 
behave statistically good. As said in Section IV, the me- 
dian of G gives a good estimation of the central value of 
its statistical distribution if any monotonic function of G 
behaves normally. Although the statistical distribution 
of the conductance becomes almost Gaussian for smaller 
degrees of disorder, I have systematically used the me- 
dian in the rest of this Section in order to be completely 
self-consistent. 

Figures 2 and 3 show the evolution of the scaling be- 
havior of the conductance as disorder diminishes. While 
finite size effects arc still unimportant for W — 5, a region 
in which conductance increases linearly with the system 
size has clearly developed for W = 2. This region ex- 
tends more and more as W approaches zero. Actually, 
the inference of a localization length beyond W = 2 is 
not possible analyzing conductance data for sample sizes 
below L = 100. As commented before, only systems hav- 
ing a localization length below £ ~ 10 4 show numerical 
indications of localization for sample sizes of the order 
of L = 100. In other words, when £ ^> L one just ob- 
tains the linear increase of conductance due to the linear 
increase of transport channels. 

Fig. 4 summarizes the main results of this subsection: 
localization length is plotted as a function of disorder. It 
seems clear that £ diverges for vanishing disorder. In fact, 
a x 2 fit by power law gives an exponent very close to —4 
and even the simpler fit by 1 + aW~ 4 works pretty good 
(see Fig. 4). Localization length values are also com- 
pared with results obtained by finite size scaling-based on 
data collected for strips of increasing width.E3 ~E3 Surpris- 
ingly, both methods do not give identical results. Present 
values of £ are somewhat larger for large disorder but 
their evolution for smaller disorder values is not so steep. 
For small disorder values, the scaling of the conductance 
of square clusters attached to semiinfinite wires yields 
considerably smaller localization lengths. I have no ex- 
planation for this discrepancy. It could be that numerical 
extrapolation of £ from finite strips has a larger uncer- 
tainty or even that the whole procedure is dominated by 
the more localized states. Let me repeat here that G is 
a statistical variable and that the scaling behavior of a 
single magnitude used to characterize the whole distri- 
bution of values can certainly depend on the choice (see 
Fig. 1 for an example of an extreme discrepancy). 

B. Random Magnetic Flux 

This part presents the main results of the paper. Let 
me begin with a detailed comparison between statistics 
based on the mean of the conductance distribution and 
statistics based on its median. Fig. 5 compares in a 
logarithmic scale the behavior of both magnitudes at an 
energy E = —3.3 close to the band limit. The values 
shown for L — 100 are not so precise as the rest because 
a sensible percentage of the samples have conductances 
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that are smaller than numerical precision. Apart from 
minor size effects at small cluster sizes, both magnitudes 
show exponentially decaying behavior. While localiza- 
tion does not offer doubts, the exact value of localization 
length depends on the magnitude used to characterize 
the central value of the distribution of conductances. £ 
increases from 20.9 for the median to 60.6 for the mean. 
As said previously, the first value is the correct one since 
the distribution of the conductance logarithm is indeed 
normal and gives the same median. Fig. 6 shows that the 
problem becomes softer as soon as energies well within 
the band are studied. Although conductance is still be- 
low one quantum unit, its distribution comes closer to a 
Gaussian distribution and therefore, mean and median 
do not differ as much as in the previous example. 

Let us analyze now the trend followed by the size de- 
pendence of conductance when energy is varied. While 
localization is clear near the band edges, it is quite dif- 
ficult to get some insight into the scaling behavior of 
conductance at energies \E\ < 3. Fig. 7 shows my 
larger data collection corresponding to an arbitrary en- 
ergy E — 0.1751 close to the band center. 5500 measure- 
ments of the conductance have been necessary to reduce 
the error bar of the mean to the small values shown in the 
Figure for L < 70. At larger sizes, a still very large num- 
ber of 500 conductances of randomly generated samples 
have been calculated. Error bars are considerably larger 
but not so much to preclude a precise analysis. In my 
opinion, the decrease of the conductance for increasing 
system sizes does not offer any doubt. Moreover, a spe- 
cial calculation at L = 200 has been carried out using a 
different algorithm that allows the study oi-iarger systems 
paying the price of longer computer times £j The decrease 
of the conductance below its value at L ~ 100 is quite 
clear: G(200) = 1.330 ±0.014 < G(100) = 1.381 ±0.005. 

The good quality of the database for this energy gives 
the opportunity of doing a reliable fit by some chosen law. 
Using Eq.(|Tl|) with four parameters (a, b, c, £), a localiza- 
tion length of £ = 5152.8 is obtained. Fig. 7 shows the 
fitted model as a thick continuous gray curve. The rela- 
tively large value of the localization length is just below 
the threshold of visibility — given the present computer 
facilities. Let me insist once again, that the data shown 
in Fig. 7 do not prove exponential decay of conductance 
although they are nicely described by an exponential law 
including finite size corrections. On the other hand, it 
is obvious that any law with properly selected parame- 
ters is able to reproduce the roughly linear decrease of 
conductance for L > 50. Nevertheless, reproducing the 
whole scaling behavior is not so trivial. For example, if 
all conductance values except the one corresponding to 
L = 200 are fitted by Eq.(|ll|) and by a rational function 
model like: 



G(L) 



aL + bL 2 



uous line of Fig. 7) and adjusting (a, b, c, d, e) (thick 
dashed line of Fig. 7) — whereas it remains almost the 
same when the exponentially decaying model is employed 
(thick continuous gray line of the same figure) . 

Actually, a large database has been also collected for 
the conductance distribution at the band center. Mean 
conductance steadily increases up to the largest calcu- 
lated size (L = 200). A very large localization length 
£ = 16707.5 is obtained using Eq.(ll) with four param- 
eters, although numerical uncertainty is as large as its 
value. The steady increase of the conductance at the 
band center of this model w.as firstly obtained in the sem- 
inal work of Lee and Fisher t3 Although they did, in prin- 
ciple, the same kind of computation, numerical values do 
not coincide for some unknown reason (the inset of Fig. 2 
of the referred paper shows conductances that are about 
seven times smaller than the ones obtained in this work) . 
In any case, either £ is infinite just at the band center or 
larger than any other calculated localization length. 

Let me point, that the mean conductance has been 
systematically used to get the localization length of this 
model for several energies within the band. Numerically, 
the advantage of the mean relative to the median is its 
better normal distribution for finite samples £3 Its use 
is justified as soon as conductance distribution behaves 
normally (Gaussian distribution). Fig. 8 shows that this 
is indeed the case at E = 0.1751 both for L — 64 and for 
L = 200. The small decrease of the mean conductance is 
not distinguishable in the distribution. This figure also 
shows the typical width of about one quantum unit that 
has been observed for the conductance fluctuations in 
numerous occasions (the standard deviation of the dis- 
tribution obtained for L = 200 equals 0.314). 

Fig. 9 collects the localization length obtained at sev- 
eral energies. The almost exponential increase of £ near 
the band edges perfectly reproduces the findings of previ- 
ous studies!] On the other hand, the localization length 
seems to reach a saturation value below 10 in the larger 
part of the band being the band center the only possible 
exception. Two comments are opportune at this point. 
First, a continuous logarithmic increase of the localiza- 
tion length up to the band center implies a more than 
macroscopic value of £ for almost all states in the band. If 
localization is much larger than the sample size, the own 
concept of localization looses its meaning. In my opin- 
ion, the scenario shown in Fig. 9 is much more attractive: 
localization length is very large but not infinite, that is, 
irrelevant. Second, although renormalization group the- 
ories also predict £ ~ Zexp(^4i 2 ) where I is the mean free 
path and A some finite constant, my numerical results do 
not necessarily contradict this behavior because the RMF 
model gives values of|-L,of order of 1 when non-diagonal 
disorder is maximum E3 



l + cL + dL 2 + eL 3 



(13) 



the fit fails to reproduce the last point if the rational 
function is used — both adjusting (a, c, d) (thin contin- 
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C. Random Hopping Sign 

The last case studied in the paper corresponds to a 
model of real hopping with a random sign. The Hamilto- 
nian is given by (Q) when link fields are limited to values 
and h. The band extends from ~ — 2\/3 to ~ 2\/3 in this 
case. Being the amount of non-diagonal disorder fixed, 
conductance can be studied as a function of energy. Since 
this model belongs to a universality class that is different 
from the previous one (orthogonal instead of unitary), 
some differences in the scaling behavior of conductance 
could, in principle, be expected. Nevertheless, numerical 
results show only minor changes. 

Scaling of conductance is shown in Fig. 10 for two dif- 
ferent energies: E = (band center) and an arbitrary 
energy E = 0.1751 close to it. While localization length 
is almost infinite (£ = 23317) in the first case, it is not 
very large in the second one (£ = 461.7). As in the com- 
plex model, an infinite value of £ is compatible with con- 
ductance data at the band center but not with the clear 
decrease shown by the conductance at the second energy. 
In both cases, the mean and its error bar have been used 
because the distribution of conductance values follows a 
Gaussian law in spite of its small average value (Mean 
conductance is of the order of 1 quantum unit when the 
conductance of an ordered sample is as large as 100 at the 
band center of a cluster of side L = 100). Conductance 
scaling at other energies follows a trend that is quite sim- 
ilar to the one followed by the previously analyzed RMF 
model. It can be concluded that the universality class of 
the model does not matter noticeable in the scaling be- 
havior of conductance although numerical values of the 
localization length are quite different (£ = 461.7 (5152.8) 
for RHS (RMF) model at E = 0.1751). 

VI. DISCUSSION 

Conductance sampling as a function of the system 
size has been methodically used to infer the localization 
lengths of a variety of disorder models. Results for the 
Anderson model of diagonal disorder are globally consis- 
tent with the ones obtained using well known approaches. 
Nevertheless, numerical values of £ are somewhat surpris- 
ing: although localization length diverges for vanishing 
amount of disorder, the divergence is weaker than the 
one provided by other methods. The origin of the dis- 
crepancy is not known although it may be due to the 
non-normal distribution of conductances when the local- 
ization regime is approached. 

Both the real — belonging to the Gaussian Orthogo- 
nal Ensemble universality class — and the complex (GUE 
universality class) models of non-diagonal disorder stud- 
ied in this work show decaying values of the central con- 
ductance as a function of cluster side at all studied en- 
ergies. When fitted with an exponential law corrected 
for finite-size effects, values of the localization lengths 



are obtained. Although large they are typically below 
10 4 , a value that is the maximum localization length giv- 
ing a conductance that diminishes within the size val- 
ues presently accessible. The behavior of £ as a function 
of energy is remarkable. After a region near the band 
edges in which £ increases exponentially, a region of al- 
most constant localization length follows. Nevertheless, 
isolated energies giving diverging values of £ cannot be 
discarded by this procedure. Actually, the band center is 
a candidate for showing a diverging localization length. 

In my opinion, this work should motivate experimen- 
tal research in the mesoscopic transport field directed to 
the determination of the localization length of disordered 
samples. As shown by the large amount of data presented 
in this work, the decay of the conductance starts at sys- 
tem sizes well below the localization length. This fact 
makes the experimental observation of localization rel- 
atively easy. Only after having experimentally related 
values of disorder (or what is equivalent, values of the 
mean free path) with values of the localization length, 
could localization theory make a major step from conjec- 
ture to fact. On the other hand, the correct microscopic 
description of the conductance scaling at a finite tem- 
perature would greatly expedite a detailed comparison 
between experiment and numerics. 
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FIG. 1. Scaling of conductance for the Anderson model of 
diagonal disorder with W = 10. The median of the conduc- 
tance distribution is shown by circles in a logarithmic scale 
while the means obtained from the same samples are shown 
by squares. Simple linear regression provides a localization 
length of f = 6.1 for the median and £ = 11.1 for the mean 
(continuous lines). 
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FIG. 2. Scaling of conductance for the Anderson model of 
diagonal disorder with W = 5. The median of the conduc- 
tance distribution is shown by circles. Fitting by Eq.([l2|) gives 
a localization length of £ = 98.8 (continuous line). 
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FIG. 3. Scaling of conductance for the Anderson model of 
diagonal disorder with W — 2. The median of the conduc- 
tance distribution is shown by circles. Fitting by Eq.([l2]) gives 
a localization length of £ = 2382.0 (continuous line). 



FIG. 4. Localization length for the Anderson model of di- 
agonal disorder. Triangles give my estimates, the dashed line 
shows the approximate law written in the box, and the contiBp 
uous line plots 1.1 exp[(ll/W / ) 2 ], the behavior given in Ref.E3. 
Localization length is measured in units of the lattice constant 
and the width W of the distribution in units of the hopping 
energy t. 
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FIG. 5. Scaling of the conductance for the Random Mag- 
netic Flux model at an energy E = —3.3 close to the band 
edge. Circles give the median of the conductance distribu- 
tion evaluated for samples including 5000 data while squares 
show the mean of the same distributions. Localization length 
changes from £ = 20.9 for the median to £ = 60.6 for the 
mean. 
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FIG. 6. Scaling of the conductance for the Random Mag- 
netic Flux model at an energy E = —3.1 not so close to the 
band edge. Circles show median values with a rough error 
estimate while squares show the evolution of the conductance 
mean. In the last case, the error bar is of the order of the 
symbol. The clear decay corresponds to a localization length 
of £ = 378.4 for the median and £ = 592.5 for the mean. 



10 




FIG. 7. Scaling of the conductance for the Random Mag- 
netic Flux model at energy E — 0.1751 close to the band 
center. 5500 "measurements" of the conductance have been 
possible up to L = 70. From this size on, only 500 define the 
average value with a larger error bar. The average conduc- 
tance at L = 200 has been calculated using a different nu- 
merical inversion subroutine. Finite size scaling with Eq. ([tl|) 
gives £ = 5152.8. 



FIG. 8. Conductance distribution of the Random Magnetic 
Flux model at energy E = 0.1751 for two different sizes: i) 
L — 64 (dotted line) and ii) L = 200 (thick continuous line). 
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FIG. 10. Scaling of the average of conductance for the Ran- 
dom Hopping Sign model at the band center (right part) and 
at an arbitrary energy E — 0.1751 close to it (left part). 
The error bars indicate the precision of the mean of the con- 
ductance distribution. The exponential decay corresponds to 
a localization length of f = 23317 at the band center and 
£ = 461.7 close to it. 



FIG. 9. Localization length for the Random Magnetic Flux 
model as a function of energy. The band limits are approxi- 
mately at ±2v3- 
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